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Summary. — The effects of magnetism on high pressure properties of transition 
metals and transition metal compounds can be quite important. In the case of Fe, 
magnetism is responsible for stability of the body-centered cubic (bcc) phase at 
ambient conditions, and the large thermal expansivity in face-centered cubic (fee) 
iron, and also has large effects on the equation of state and elasticity of hexagonal 
close-packed (hep) iron. In transition metal oxides, local magnetic moments are 
responsible for their insulating behavior. LDA-I-U results are presented for CoO 
and FeO, and predictions are made for high pressure metallization. The inclusion 
of a local Coulomb repulsion, U, greatly inhibits the high-spin low-spin transitions 
found with conventional exchange-correlation functionals (i.e. generalized gradient 
corrections, GGA). We discuss theory and computations for the effects of magnetism 
on high pressure cohesive properties. 



1. Magnetism 

One might think that one needs to be concerned with magnetism only if one is inter- 
ested in magnetic properties of materials. This is not the case! Rather the stable crystal 
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structure can depend on properly taking into account magnetism, which can strongly 
affect phase stability, lattice distortions, elasticity, equations of state, and vibrational 
frequencies. Furthermore, magnetism can be important in materials containing tran- 
sition metal ions even when the magnetic moments are not ordered. The distinction 
between "non-magnetic" materials and disordered, paramagnetic, materials is crucial. 

Magnetism in crystals comes about because electrons have a quantity called spin, 
which is a vector quantity that behaves like quantum angular momentum. Electrons 
also have orbital angular momentum, which can lead to orbital magnetic moments. In 
heavy transition metal ions with partially occupied f-electron shells, such as rare earths 
or lanthanides, these orbital moments can be very important. 

Here we will concentrate on d-electron systems, but there is much complex and fas- 
cinating physics in the f-electron systems. In the f-electron systems, the f-electrons tend 
to form very localized, strong magnetic moments on each atom. The localized f-electrons 
then interact with delocalized band- like states, and this interaction can lead to interest- 
ing phase transitions with pressure, for example in Ce [1]. Electronic d-states tend to be 
more delocalized than f-states, though localization is important, as we will discuss below. 

The fact that electronic wave functions are antisymmetric with respect to exchanging 
two electrons leads to the Paiili exclusion principle, without which electrons would all 
fall into the nucleus, and there would be no atoms. The exclusion principle states that 
two electrons cannot be in the same state. If electrons did not have spin the universe 
would be a very different place, since only one electron could then be put into each band. 
Since electrons do have spin, it is possible for two electrons to occupy the same state, 
with one electron having the opposite spin of the other. In many atoms, molecules, 
crystals, and liquids, the electrons are paired up, with each member of the pair having 
the same probability distribution in space, but with opposite spins. These spin states 
are conventionally called "spin-up" and "spin-down", though there is nothing "up" or 
"down" about them except the way they might be drawn. Note that one must be careful 
to Tinderstand that the concept of electron pairing in states is a rather loose way of 
talking. In reality electrons are indistinguishable from each other, and one should talk 
about the quantum states being paired, rather than particular electrons. 

Electrons have an effective interaction that is different depending on whether they 
have the same spin or not. Since electrons with opposite spin can occupy the same space, 
they have a higher potential energy of interaction on average. Thus all else being equal, 
electrons would want to have the same spin in order to lower the system's potential 
energy. A system with electrons that have the same spin direction is a ferromagnet. 
There is no free lunch, though. The cost of lowering the potential energy by lining up 
the spins is to raise the system's kinetic energy, since higher states must be occupied 
instead of the lower energy paired states. Thus there is competition between electronic 
potential energy which favors magnetism, and the electronic kinetic energy, which favors 
a non-magnetic electronic structure. As pressure is increased, electrons are pushed close 
together, and the relative potential energy change between paired and unpaired electrons 
becomes less important; bands become wider, making the kinetic energy cost smaller, 
so that in general materials become non-magnetic with increasing pressure. The total 
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energy change between the magnetic and non-magnetic system in the simple picture 
presented so far is known as the exchange energy. 

As temperature is raised from low temperatures, the magnetic-moment directions on 
each atom will fluctuate more and more, and at some critical temperature, called the 
Curie temperature, or Tc in ferromagnets, the moments will disorder. In general, there 
are still magnetic moments on the ions above Tc, they are just disordered in direction. 

Antiferromagnets have moments of opposite direction on alternating sites. It is the 
different hybridization of electronic states that leads to antiferromagnetic rather than 
ferromagnetic order, so that the kinetic energy is lowered. This is sensitive to pressure, 
so some ferromagnets become antiferromagnetic with increasing pressure, as in fee iron. 

In some cases the lowest free energy state has non-collincar spins. This can arise from 
"frustration," which is the situation where it is impossible to tile a lattice with a per- 
fect antiferromagnetic pattern, with every atom having only neighbors with the opposite 
pointing spin. Examples are fee and hep structures. Such effects can be very impor- 
tant to material properties, and are responsible for the anti-Invar effect (high thermal 
expansivity) in fee Fe, for example. 

This is not meant as a comprehensive review, but rather as lecture notes, and an 
introduction of this complex field to the student. Examples are mainly chosen from our 
own theoretical work, or from published experiments. 

IT. Itinerant magnetism. - There are two endmember models for understanding an- 
tiferromagnetism. In the band, or Slater, picture [2], it is the different exchange inter- 
actions between like- and unlike-spin electrons, combined with the kinetic energy of the 
resulting different band states that leads to the stability of antiferromagnetic, rather than 
ferromagnetic or non-magnetic states. Antiferromagnetism is a zone-boundary instabil- 
ity that leads to a doubling of the unit cell, giving folding of the electronic states. These 
states then hybridize differently through the exchange potential than they would have in 
the ferromagnetic case. 

The Stoncr model demonstrates the band picture of magnetism, and though a simple 
model, turns out to be predictive. In the Stoner model, the magnetization energy AE is 

-/M2 M2 

where M is the magnetic moment or magnetization, the Stoner Integral / is an atomic 

property, A^(0) is the density of states at the Fermi level (or top of the valence band). The 
first term on the right hand side is the exchange energy and the second is the change in 
the band energy with magnetic moment. Minimization of AE gives the Stoner criterion: 

(2) IN{0) > 1 



for a magnetic state to be stable. In the more sophisticated extended Stoner model, the 
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average density of states, 



(3) 



N{M) = M/Se, 



is used, where Se is the spin (exchange) splitting, and 



(4) 



iV(0) = N{0) = dM/d6e. 



The instabihty criterion is IN{M) > 1. As pressure is increased the effective density of 
states generally decreases, and magnetism will decrease and disappear with increasing 
pressure. 

Figure 1 shows Stoner diagrams for CoO, hypothetical FeSiOa in the cubic perovskite 
structure, and FeO in the NiAs structure. With increasing pressure the effective density 
of states drops, and the stable magnetic structure becomes low-spin due to this effect. 

1"2. Mott insulators. - Predicting properties of transition metal-bearing oxides is 
a severe problem in modern band theory, although qualitative understanding of the 
behavior of transition metal oxides has been developed [3, 4]. Examples of problem 
compounds are CoO and FeO; conventional band theory (the local density approximation, 
LDA, or generalized gradient approximation, GGA [5]) predicts them to be metals, but 
they arc actually insulators. Materials which are insulating because they have magnetic 
correlations are known as Mott insulators [6]. In contrast to density functional theory 
(DFT) calculations within LDA or GGA, Hartree-Fock theory gives a large band gap 
for transition metals oxides [7], much larger than the experimentally observed gaps. It 
is generally understood that the the problem with LDA and similar theories for these 
materials is that they underestimate the local Coulomb repulsions between electrons. 
Band theory assumes that electrons are delocalized through space, but if the electrons 
can be considered as localized, the energy will depend on how many electrons are on a 
given site (i.e. localized on a given atom or region of space). In order to conduct electricity 
the energy cost for an electron to move from site to site cannot be too high. There are 
several approaches that give a insulating behavior for the problematic transition metal 
oxides. 

In the simplest one-band Hubbard model, the Hamiltonian is given by 



and the band states are split into an upper and lower Hubbard band, with a splitting of 
U. The first term is the hopping, or hybridization, governed by T. The creation operator 
c, adds an electron to site i, and the destruction operator cj subtracts an electron from 
site i. The number of up electrons on site i is 



(5) 




(6) 
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Fig. 1. - Stoner diagram for (a) CoO in the NaCl (Bl) structure with a cubic lattice, (b) FeSiOs 
in the cubic perovskite structure, and (c) FeO in the NiAs structure. The effective density of 
states is plotted against moment for several different volumes. The inverse Stoner parameter, 1/J 
is also shown, and places where the effective density of states crosses 1/7 with a negative slope 
arc stable solutions. At high volumes (low pressures) there is only a high-spin (large moment) 
solution. At intermediate pressures both high-spin and low-spin solutions exist. (Whether such 
coexistence is possible, i.e. whether the transition is first-order or second-order, depends on the 
shape of the effective density of states.) At high pressures, there is only a low-spin solution. The 
Stoner model behaves very differently from the localized, crystal field, picture (shown in inset) 
which predicts a discontinuous transition, never giving intermediate moments. Prom ref. [8]. 
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and similarly for the down electrons. In this simple model there can be 0, 1 (up or 
down), or 2 (one up and one down) electrons on a site. Even this simple model has never 

been solved exactly in three dimensions. The physical picture, however, is clear. The 
first term is the band term. The second is a local repulsion. A local repulsion between 
electrons can open up a gap in the excitation spectrum, and make a insulator out of a 
band metal. 

1'2.1. LDA+U. One approach that gives an insulating ground state for transi- 
tion metal oxides is the LDA+U model, which adds a Hartree-Fock-like local Coulomb 
repulsion tensor U and exchange interaction J, and then attempts to correct for dou- 
ble counting. The LDA+U model has given excellent results for a variety of systems 
[9, 10, 11, 12], but the limits of the model are still not well known, for example whether 
it properly predicts the high pressure behavior of Mott insulators. Modern LDA+U is 
rotationally invariant, but still has some dependence of the choice of local orbitals in 
which to apply the corrections. The rotationally invariant LDA+U contribution to the 
energy is given by 



(7) 



2 

mim'^m2m2<T 



where Edc is the double counting correction accounting for the on-site Coulomb inter- 
action already included in LDA. The tensors U and J are the Coulomb and exchange 
integrals between electrons in orbitals m^. The n's are the site occupancy matrices and 
a designates spin. For the double counting correction, we use the form: 

(8) Eff"'^^^ = ^U2n{2n - 1) - ^J[n^n^ - 1) + n^{n^ - 1)]. 



where 



^ " (21 + IV ^ ("^"^'1^1™™') 



^ 2/(2/ + 1) 



; 7 'S^ ((mm'\ — \mm') — (mm'l — lm'm)) 



(9) 



and where n'^ = J2m "'^mm^ ^^'^ ^ = (^^ + n^)/^- The Coulomb and exchange tensors U 
and J are defined by 

Umim2m{m'2 = (miwi | ^ Imama) 
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(10) Jm^m^m'^m'^ = (mim'i | ^ |m2m2), 

and are to be evaluated over localized orbitals. In practice U and J are input parameters, 

but can be determined from constrained occupancy computations [12, 13]. Actually the 
Slater integrals Fq, F2, and F3 are input for d-states, and U = Fq, J = F2 + F^/IA, and 
F2/Fi = 0.625. We use the atomic values for F2 and F3: ^2=8.18 eV and F4=5.15 eV, 
or J=0.95 eV, for CoO and F2=7.67 eV and ^4=4.79 eV, or J=0.89 eV for FeO. 

LDA+U has been shown to give good predictions of the electronic structure of NiO 
[12, 14, 15]. In ref. [14] the equation of state and elastic constants for cubic NiO were 
computed with LDA+U and SIC (see below), and reasonable agreement with experiment 
was found. However, there have been few tests of LDA+U total energies, and here we 
study in detail the energetics with respect to strain and compression in FeO. 

1'2.2. Self-interaction corrections. Another successful model, which still has not 
been fully explored for Mott insulators, is known as the self- interaction correction (SIC). 
This is understood easiest by considering a hydrogen atom with one electron. The self- 
consistent field electrostatic energy and Hartree self-consistent potential are computed 
from the electronic charge density. The charge density for a hydrogen atom is spread 
out, yet actually an electron is only at one place at a time; there should be no electron- 
electron interaction when there is only one electron. The Hartree potential V^and energy 
are computed from the charge density, so that 



(11) V„(r) = j 



for example, which is an electron-electron interaction, which should not be present in 

a hydrogen atom! The LDA exchange-correlation potential also depends on the charge 
density, and would be non-zero in the hydrogen atom. In the exact density functional, 
which is unknown, the exchange-correlation potential must exactly cancel the Hartree 
potential, giving zero for the total electron-electron interaction. This can be enforced, 
by correcting the potential so that 

(12) Vsic = Vlda -Y^I dr^V {piu) 

where 



(13) Pik = iJtkAk- 

The simple SIC correction above clearly depends on the localization of the basis 
set. It works very well in an atom, but in a crystal the Bloch functions are extended 

throughout space, and the correction goes to zero. However, one can perform a unitary 
transformation on the orbitals and find sets of localized orbitals. It has been shown that 
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there is a variational principle, and the most localized orbitals have the lowest SIC energy. 
Thus SIC computations are quite computationally intensive, as there is an inner iterative 
loop in which the orbitals are localized. SIC does very well for excitation energies and 
does predict an insulating ground state for the transition metal oxides [16, 17, 18], and 
gives a qualitatively correct picture for the rare earth metals [19]. 

In the exact density functional theory, there should be no self-interaction. One can 
still question whether the way SIC enforces freedom from self-interaction is realistic or 
accurate. For example, the weighted density approximation (WDA) is constructed to 
give no self-interaction per orbital, yet it makes only small changes in the electronic 
structure in most systems studied so far [20]. In contrast, SIC makes large changes in 
electronic structure for all systems studied. 

LDA-f U and SIC make very different predictions for the electronic structure. Whereas 
LDA-I-U pushes unoccupied states up in energy, SIC pulls occupied states down in energy. 
This makes a difference when there are different types of bands present. For example, 
in the transition metal oxides, 4s states are not affected directly by U. The nature of 
the band gap is thus different in the two models, as will be discussed further below. 
In ref. [14] the EELS spectra are computed for O 2p for LDA-I-U and SIC, and both 
compare well with each other and with experiment. However, the differences in models 
are minimized by considering only the O 2p states, which are only indirectly affected by 
the correlations. 

1'2.3. Dynamical mean field theory. Dynamical mean-field theory (DMFT) is a 
method for treating correlated systems that becomes exact in infinite dimensions [21]. 
It can be applied as an extension of LDA+U that includes a frequency dependent hy- 
bridization function [22]. A multiplet ion is solved self-consistently within the mean 
field Anderson impurity model. It is very time-consuming with few cases studied, and 
it is probably not currently tractable to study a system such as FeO using DMFT as a 
function of strain and pressure, as we did here. In DMFT, one still uses the parameter 
U. 



1"2.4. Hartree-Fock. In contrast to density- functional-based methods, Hartree-Fock 

theory gives a large band gap (way too large) for transition metal oxides, and some 
feel that Hartree-Fock should therefore be taken as the zeroth order method for these 
materials [4, 23, 7]. This is indeed the basis of the LDA-I-U method. Hartree-Fock by 
itself is not a reasonable way to study the high-pressure properties of Mott insulators, 
because it would grossly overestimate the pressure of a metal-insulator transition, since it 
greatly overestimates the gap. This is because the Couloumb repulsion in Haretree-Fock 
is completely unscreened. Furthermore, Hartree-Fock does not work well for any metal, 
always predicting a singularity in the density of states at the Fermi level. Thus it is not 
a good way to study metal/insulator transitions. 
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2. - Methods 

We have used a variety of methods to study Fe, FeO and CoO, and only a brief 
outhne of the methods will be given here. Computational details are given in the ref- 
erenced papers. For Fe our most accurate method is the Linearized Augmented Plane 
Wave (LAPW) method. We have used the LAPW method to study non-magnetic and 
collinear magnetic properties of Fc [24, 25, 26]. LAPW is a full-potential all-electron (i.e. 
no pseudopotential) method, and the basis is very flexible, and is suited both to the inter- 
stitial region and the atomic cores. One advantage of LAPW is that it is straightforward 
to converge the results with respect to basis set size. Our computations arc done using 
the Generalized Gradient Approximation (GGA) for the exchange-correlation potential 
[5]. 



3. — Results and Discussion 

3'1. Overview of effects of pressure on magnetism. - A straightforward effect of pres- 
sure on magnetism is through structural phase transitions. Compression can drive struc- 
tural phase transitions, and the magnetic properties of the different phases will be differ- 
ent. For example, body-centered cubic iron is ferromagnetic, and hexagonal close-packed 
iron was believed to be non-magnetic, so the bcc to hep transition would also be a ferro- 
magnetic to non-magnetic transition. One can also have transitions from one magnetic 
structure to another. The simplest is the Curie point, where the magnetic moments 
become disordered. Generally, there are still local moments on the atoms, but they are 
no longer ordered. The Curie temperature in bcc iron is 1043 K. One could also have a 
transition from ferromagnetic to antiferromagnetic, or from an ordered structure to an 
incommensurate structure. 

As pressure is increased magnetic moments tend to decrease, and eventually mag- 
netism is squeezed out. This can be understood in the Stoncr model as due to the 
general increase in bandwidths with pressure, decreasing the effective density of states, 
whereas the effective Stoner parameter is approximately constant. This also happens in 
the Hubbard picture, since U is approximately constant or decreasing with pressure, and 
the hopping parameter T increases in magnitude with pressure. Actually the parameters 
of the two models are closely related. Actual computations are shown for FeO using the 
LAPW method and the GGA in fig. 2. Note that the magnetic behavior with pressure 
depends strongly on the crystal structure (cubic lattice vs. one with a small rhombo- 
hedral strain) and on the magnetic order (antiferromagnetic vs. ferromagnetic). The 
antiferromagnetic state with cubic lattice shows a strong first-order high-spin low-spin 
transition with an appreciable AV of 8%. With a rhombohedral strain, or with ferromag- 
netic order, the transition is gradual. This kind of comparison is one of the useful things 
that can be done with theory. In experiments one is generally stuck with the lowest free 
energy phase unless there is a large activation energy so that metastable phases can be 
studied. With theory one can answer the question. "What if FeO were ferromagnetic?" 
High-spin low-spin transitions are predicted in all the transition metal oxides [8, 27], and 
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Fig. 2. - Computed magnetic moments for FeO using the LAPW method and the GGA. For the 
cubic lattice the antiferromagnetic solution gives a discontinuous high-spin low-spin transition. 
However, with the equilibrium rhombohedral distortion at each volume, the moments decrease 
continuously. The behavior also depends on the nature of magnetic ordering, i.e. antiferromag- 
netic or ferromagnetic. 



have been observed in sulfides and f- metals [28]. 

Metals such as Co and hep Fe do not show discontinuous transitions in our compu- 
tations, but show a gradual reduction in moments as the bands widen with increasing 
pressure (fig. 3) [26]. The eflFect of pressure on magnetism is a strong function of the 
magnetic structure, as well as chemical composition. 

3'2. Magnetic behavior of Fe, FeO, and CoO with increasing pressure. - 

3'2.1. Fe. There are three well-known crystalline phases for Fe: body-centered cubic 
(bcc, or a-Fe), which is the stable form at ambient conditions, face-centered cubic (fee, 
or 7-Fe) which is stable at high temperatures, and hexagonal close-packed (hep, or e- 
Fe), stable at high pressures (fig. 4). There is also a small bcc field ((5-Fe) just before 
melting at low pressures. Bcc iron is ferromagnetic, fee has magnetic correlations, but 
not an ordered magnetic structure, and hep was believed to be non-magnetic (though, 
see discussion below). See ref. [29] for a review of experimental studies of Fe at high 
pressures. 
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Fig. 3. - Computed magnetic moments for afml (open symbols) and afmll hep Fe. These are 
the moments in the muffin tins, which have a radius of 2.0 bohr. Note that in the original 
presentation of these results (ref. [26]) the moments were plotted incorrectly smaller by a factor 
of two. 



One of the well-known failures of the local density approximation (LDA) is that it 

gives the wrong ground state for Fe. LDA only gives the magnetic bcc phase at negative 
pressures, but at zero pressure gives a close-packed non-magnetic ground state. One of 
the early successes of the GGA was the correct prediction of the bcc ground state and 
an accurate transition pressure to hep from bcc of 11 GPa [24, 30] compared with a 
10-15 GPa hysteresis loop from experiment. This showed that the GGA was accurate 
for the magnetic stabilization energy of bcc-Fe. Furthermore, bcc iron is only stable due 
to its magnetism. Non-magnetic bcc iron would not be a stable phase, except possibly 
in the S-Fc field just before melting. There had been speculations that the Earth's solid 
inner core was bcc-iron, but calculations showed that bcc-Fe was mechanically unstable 
at those conditions because it was non-magnetic at such high pressures [25]. 

Not only is magnetism important in bcc iron, but it is also important in fee. We tried 
to find the phase transition from hep to non-magnetic fee using the particlc-in-a-cell 
model (ref. [31] and unpublished), but did not find a stable field for fee below melting. 
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Fig. 4. - The experimental phase diagram for Fe. See ref. [29] for discussion. 



which could be explained by magnetic stabilization of fee. Experimentally, it is known 
that fee iron is magnetic (i.e. has local moments). 

The computed equations of state for magnetic bcc-Fe, and non-magnetic fee and 
hep, are shown in fig. 5. The inset shows the total energies for these three phases; the 
transition from bcc to hep occurs at the common tangent. The computed P-V equation 
of state is in good agreement at high pressures, and with bcc, but discrepancies are seen 
at lower pressures; the non-magnetic hep equation of state is too stiff at low pressures. 
The discrepancy at low pressures seems very large if one compares the zero-pressure bulk 
modulus Kq from the equation of state. The experimental value is 165 GPa, compared 
with 292 GPa from the GGA equation of state. Kq is a Active quantity for hcp-Fe, 
since it has not proved possible to quench hcp-Fe to zero pressure, but nevertheless 
the discrepancy seems larger than expected compared with GGA results for other hep 
transition metals (for example for Co Ko=190 GPa from experiment and 212 GPa from 
GGA, and Re is 365 GPa from experiment and 344 GPa from GGA [26]). 

Disagreements between theory and experiment often lead to advancement in theory, 
or in our understanding of physics. In this case, it appears that the discrepancy in the 
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Fig. 5. - Computed equation of state for Fe using the LAPW method and GGA. A transition 
pressure from bcc to hep of 11 GPa is obtained from the common tangent of the energy curves 
shown in the inset, in good agreement with experiment. The computed PV equation of state is 
in good agreement with experiments [32, 33] at high pressures, and with bcc. 



equation of state of hcp-Fe is not due to inaccuracy of the theoretical calculations or 

methods, but is due to the neglect of magnetism in hep-Fe. Experimental data had 
been interpreted to show no magnetism in hcp-Fe, and indeed calculations showed that 
ferromagnetism was not stable in hcp-Fe. Mossbauer experiments showed no evidence 
for ordered magnetism in hcp-Fe, even down to 0.03 K [34, 35]. However, even in the 
earliest papers on Mossbauer in hcp-Fe it was recognized that magnetism was not ruled 
out by the data[36, 37]. Given the large unexplained discrepancy in the equation of 
state, we started looking for stable magnetic structures. Two stable antiferromagnetic 
hep structures were found, and the most stable, afmll, is stable up to about 50 GPa. The 
computed local magnetic moments are shown in fig. 3. The afml structure consists of 
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layers of Fe with alternating spin along the c-axis. The afmll structure (fig. 6) alternates 
spin-up and spin-down layers along the hexagonal a-axis [26] . 




Fig. 6. - Antiferromagnetic groundstate of hep iron (afmll). Filled symbols show the atomic 
positions at z=l/4, the open symbols at z=3/4 with the arrows indicating the direction of spin 
on the atoms. The spacegroup of the afmll structure is Pmma with the atomic positions of 
the spin up states at (1/4,0,1/3) and spin down states at (1/4,1/2,5/6). Also shown are the 
orthorhombic unit cell vectors in the x-y-plane (6 and c): we chose the orthorhombic a-axis 
along the hexagonal c-axis (out of the plane). The orthorhombic b axis then coincides with 
hexagonal a. The orthorhombic b and c also define the eigenvectors for the displacements of the 
zone center TO modes (TOb and TOc) in the afmll structure. 

The bulk modulus Kg for afmll hcp-Fe is 209 GPa, a vast improvement from the non- 
magnetic GGA value of 292 GPa, but higher than the experimental value of 165 GPa. 
Another piece of evidence for local antiferromagnetic correlations comes from Raman 
spectroscopy. The hep structure has a single Raman active mode. Raman experiments 
on hcp-Fe show a second broad peak [38] which cannot be explained within hep symmetry. 
The afmll structure has magnetic symmetry lower than hep, is orthorhombic, and has 
two transverse optic Raman modes. It is very interesting that the predicted splitting in 
Raman frequencies due to the antiferromagnetic order for afmll is in excellent agreement 
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with experiment (fig. 7). However, the second peak seen in experiments is quite broad, 
suggesting that a long-range ordered magnetic structure is not present. 
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Fig. 7. - Transverse optical frequencies as a function of atomic volume. Non-magnetic calcula- 
tions are shown in filled squares. The afmll structure has two transverse optical modes (filled 
circles) with TOb being the lower and TOc the upper branch. The dashed lines through TOb 
and TOc are finite strain fits to the results to third order in V~'^^^ . Experiments [38] identify 
two peaks in the Raman spectra up to 40 GPa. The stronger, low frequency peak is shown 
with triangles up, the weaker, high frequency peak with triangles down. The inset compares the 
calculated split in TO frequencies (circles) with the Raman experiment (triangles). 



There is increasing theoretical evidence that hep-Fe has local magnetic moments, 
which are not long-range ordered, but have strong local antiferromagnetic, non-coUinear, 
correlations. On the other hand, experiments suggest that hcp-Fe is non-magnetic. 

Mossbaucr data show no ordered moments, requiring a correlation time less than 10~^- 
10~^ seconds. X-ray absorption experiments [39] also show significant decrease of moment 
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Table I. - Comparison of experimental and theoretical values of equilibrium volume (Vo)and 
bulk modulus (Ko ) for e-Fe. 

Fe(GGA) Vb (Bohr^) Kq (GPa) 

Expt[33, 32] 75.4 165 

Non-Magnetic[26] 69.0 292 

Collinear (afmll)[26] 71.2 209 



between bcc and hep, but they do not prove complete loss of magnetism; the change in 
the absorption spectrum is due to changes in density of states as well as changes in the 
spin-rclatcd satellite. Recent experiments that show superconductivity in hcp-Fc [40] 
are fascinating. (However, the data reported in ref. [40] do not prove superconductivity 
unambiguously; they could alternatively be interpreted as a magnetic phase transition. 
But they are strong evidence for superconductivity given the strong dependence of the 
resistivity of applied magnetic field.) It used to be thought that superconductivity and 
magnetism were incompatible, but that is now known to be not always true, and weak 
magnetism may promote exotic superconductivity [41]. Nevertheless, the fact that only 
specially treated samples were superconducting suggests that the superconducivity in Fe 
may be exotic. The main problem with the idea that hep Fe is locally magnetic are the 
Mossbauer measurements. It is hard to understand how at the lowest temperatures the 
spin dynamics would still be so rapid as to show no magnetism. On the other hand, it 
is rare for band theory to incorrectly predict a magnetic ground state. Including local 
correlations neglected in GGA should further promote magnetism. If there were struc- 
tural distortions or very soft vibrational modes that invalidate the Born-Oppenheimer 
approximation there would be some grounds for considering the theoretical predictions 
less firm. But hep is a very simple, close-packed structure, so that an incorrect predic- 
tion of a magnetic ground state does not seem reasonable. Furthermore, as shown above, 
including magnetism greatly improves the equation of state of Fe (table I). So there 
remains a problem reconciling the experimental and theoretical evidence. 

3'2.2. FeO and CoO. Transition metal oxides like FeO and CoO present an even more 
challenging problem to both theory and experiment. Understanding materials such as 
FeO is one of the frontier problems in condensed matter physics. FeO and CoO are Mott 
insulators, that is they are insulating because of local magnetic moments. Conventional 
band theory makes FeO and CoO metals [42] , and no small change in exchange-correlation 
potential will make them insulators. (Large changes in the exchange correlation potential 
can make them insulators, but at the expense of accurate total energies [43].) In spite of 
this failure, conventional LDA or GGA calculations predict energetic properties, such as 
equations of state, and the magnetic moments reasonably well [44]. LDA-I-U predicts a 
band gap and the canted magnetic moments and lattice strains experimentally observed 
in CoO are reproduced [45]. 

We have performed LDA-f-U computations on CoO and FeO and get very encour- 
aging results, but only with more experiments will the predictive power of LDA-I-U be 
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understood. Since the results depend on the value of U, how U is computed or estimated 
is critical. Without comparison with experiment it is hard to test different models for 
U, so we have used several different values. A complication of LDA+U is that different 
self-consistent results can be obtained depending on electronic symmetry and d-state 
occupations. Thus one can have electronic symmetry that is lower than the lattice or 
magnetic symmetry due to the orbital ordering. We report results on the lowest energy 
states we have found, but no systematic study of the different metastable states has yet 
been made. 

The equations of state for CoO obtained using LDA+U (and the GGA, i.e. U = 0) 

are reasonable (fig. 8), and a U from 2-5 eV is consistent with the experimental equation 
of state. Figure 9 shows the computed magnetic moments for antiferromagnetic CoO 
with a cubic lattice, computed using LMTO-ASA with GGA. The curve for [/ = is the 
GGA result. Note that our LDA+U results used GGA for the density functional, so our 
results could properly be called GGA+U. GGA gives a high-spin low-spin transition, as 
we found earlier [8]. 

The band gap in CoO is predicted to initially increase with pressure, and then de- 
crease. At zero pressure, the lowest conduction states are 4s states, so that the gap is an 
intra-atomic gap between 3d and 4s. As pressure increases, the 4s states arc driven up 
and the gap becomes a d-d gap, which then decreases with increasing pressure. LDA+U 
pushes unoccupied states up into the conduction band, so that the conduction band 
edge contains 4s character, which are not affected by U. The d-s nature of the gap 
was pointed out in ref. [46]; they suggested that this is consistent with photoemission 
data. In contrast, SIC lowers the occupied states, rather than raising the unoccupied 
states, and gives a prediction of d-character for the lowest conduction band states [17]. 
This suggests that LDA+U may be a better approximation for the nature of the gap in 
transition metal oxides than SIC. 

Of critical importance for gap closure, that is metallization pressure, is the behavior 
of U with pressure. The assumption is generally made that U is about constant with 
pressure, being a local, atomic-like property, but it may decrease with strong compression. 
In that light, the U =2 eV result of 170 GPa would be a lower bound for the metallization 
pressure for cubic CoO. Lattice distortion may increase this to higher pressure, as would 
a higher value of U. It is interesting that such moderate values of U are sufRcient to 
give an insulating state over a wide range of compression, when the GGA itself gives a 
metallic band structure. 

At high pressures and temperatures, FeO transforms from the rhombohcdrally dis- 
torted rocksalt (Bl) structure to a superlattice of NiAs and anti-NiAs structure [47, 48, 
49]. However, at room temperature, FeO can be maintained in the distorted Bl structure 
to over 120 GPa [50]. Here wc discuss the behavior of FeO in the distorted rocksalt struc- 
ture, since understanding the behavior of transition metal oxides in this simple structure 
is a prerequisite to understanding the behavior of more complex phases. We have per- 
formed a large number of self-consistent computations using the full potential LMTO 
method [22] for different lattice strains and d-orbital occupancies. The lowest energy 
state we find has rhombohedral symmetry, except for very high pressures with U = 4.6 
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Fig. 8. Computed equation of state for CoO and equation of state parameters from a Vinet 
equation [51] fit to the computed energies. The experimental values for the zero pressure volume 
and bulk modulus are 261 bohr^ and 181 GPa, respectively. 



cV (fig. 11). This value of U was obtained by computing the change in eigenvalues with 
potential shift in each d-orbital, and is probably the best estimate at zero pressure [13]. 
As mentioned above, the effects of pressure on U are still unknown. Smaller values of U 
{U < 3.5 eV) give a metalHc band structure at zero pressure. 

The GGA (U = 0) density of states is shown in fig. 12 at the experimental zero 
pressure volume. The electronic structure is metallic in the GGA. With a rhombohedral 
strain, some minority spin states move down, making it possible to open a gap, but still 
a gap docs not form. LDA+U does open a gap (fig. 13), as was shown in the original 
LDA+U work (ref. [9]). LDA+U also does a much better job for an isolated Fe ion in 
MgO, indicating that it might work for arbitrary transition metal ion concentrations [52]. 

Previous work has not examined the LDA+U total energies. In order to understand 
how predictive arc the total energies computed with LDA+U, wc studied the energy as 
a function of rhombohedral strain for FeO. The lattice strain A in terms of the rhombo- 
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Fig. 9. - Computed local moments for Co in cubic CoO. U=0 is the CCA result, which shows 
a first-order high-spin low-spin transition. Note that turning on U inhibits the transition. 



hedral strain parameter 6 is given by 
(14) A = (1 + 36)-^ 



1 + S S 6 
6 1 + 6 6 
5 5 1 + 5 



We find that the lattice strain and its pressure dependence are predicted better by 
LDA-I-U than by GGA or LDA (fig. 14). This is strong evidence that LDA-I-U is making 
the right kind of corrections since the energetics of strain are quite subtle. 

We also found more than one stable state for FeO with different orbital occupancies. 
The lowest energy state at low pressures was the strained lattice with rhombohedral 
symmetry. At high pressures, when U >4.6 eV, we find a phase transition to a state 
with monoclinic electronic symmetry (fig. 11). 

The band gaps for different symmetry solutions are shown in fig. 15. The lower 
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Fig. 10. - Computed band gap for cubic CoO using LDA+U. 



energy rhombohedral symmetry d-state occupancies also have a larger band gap, and 
gap closure occurs at about 250 GPa for {7=4.6 eV. The rhombohedral strain lowers 
the gap appreciably. It is interesting to note that the monoclinic structure gap closes at 
much lower pressures. This suggests that large non-hydrostatic stresses might promote 
gap closure at lower pressures than under hydrostatic conditions. 

4. — Conclusions 

Magnetism is an important contribution to the high pressure properties of many mate- 
rials containing transition metals. For pure Fe, conventional electronic structure methods 
within the density functional theory using modern exchange-correlation functional, such 
as the generalized gradient approximation (GGA), successfully predict structural and 
magnetic properties, including equations of state, phase transition pressures, elasticity, 
etc. There are still some areas with discrepancies, such as the elastic properties of hep Fe, 
but these differences are probably due both to experimental difficulties, and to not yet 
having the proper magnetic structure in the first-principles computations. Discrepancies 
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Fig. 11. - Energy difFcrcncos between various LDA + U solutions and the unstrained rhombo- 
hedral solution (dotted line) for (a) U = 4.6 eV, (b) U = 6.0 eV, and (c) U = 8.0 eV. Solid 
line, strained rhombohedral solution; dashed line, strained monoclinic solution; dash-dot line, 
unstrained monoclinic solution. 
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Fig. 12. - Densities of states for FeO at the zero pressure volume for (a) cubic lattice and (b) 
optimal rhombohedral strain computed using GGA {U = 0). Right hashing, spin-up states in 
one Fe mufhn tin; left hashing, spin-down states. 



in theoretical computations of the melting curve [53, 54] are probably due to the difficul- 
ties of sufficient accuracy in computing the solid and liquid free energies rather than any 
fundamental problem with the underlying theory. The main remaining hurdle for Fe is 
to be able to do computations with sufficient speed and accuracy, self-consistently within 
the GGA, to compute free energies of all of the phases as functions of V, T, strain, etc., 
in order to obtain a thermal equation of state, elasticity, phase transitions, magnetic, 
and vibrational properties over a wide pressure and temperature range, though much 
progress has been made [53, 54, 55, 56, 57]. Theory suggests that magnetism in Fe is 
key to understanding the high pressure behavior of Fe, for pressures at least up to 50 
GPa. On the other hand, there are no experiments that show ambiguously the presence 
of local moments in hep Fe, and many experiments imply the opposite. If indeed iron 
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Fig. 13. - Density of states for FeO with LDA+U, and [7=4.6 eV. (a) Rhombohedrally strained 
with rhombohedral electronic symmetry at GPa; (b) rhombohedrally strained with electronic 
monochnic symmetry at GPa; (c) rhombohedrally strained with rhombohedral electronic sym- 
metry at 180 GPa; (d) rhombohedral symmetry with a cubic lattice at 180 GPa. Right hashing, 
spin-up states in one Fe muffin tin; left hashing, spin-down states. 



is non-magnetic, contrary to our best theoretical calculations, important changes to our 
best exchange-correlation functionals are required. Alternatively, perhaps defects and/or 
thermal disorder are giving rise to loss of moments, as discussed above. 

The situation for transition metal oxides is not so clear. Models are available now 
that properly give insulating behavior for FeO and other transition metal oxides, such 
as LDA-I-U and SIC. It is not yet known how predictive these methods are since the 
experimental data are not yet available, though computations agree well with present 
data. Theoretical predictions of high pressure behavior within SIC, or within other 
methods such as dynamical mean field theory are also not available. Here we presented 
predictions of behavior from LDA+U that can be tested when more experimental data 
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Fig. 14. - (a) Energy vs. rhombohedral strain S for FeO at V=137.33 bohr'^/formula unit, the 
experimental zero pressure volume. The LDA+U results are an improvement over GGA. Good 
agreement is found between full potential LAPW and full potential LMTO, and little difference 
is found between LDA and GGA at fixed volume, (b) Optimized strain for FeO as a function 
of volume. Asterisks, experiment [50]. 



become available. In any case, it is clear that local magnetism is responsible for the 
insulating behavior, and is also probably key for accurately understanding lattice strains 
and elasticity. 

This is a forefront area, and we expect to see many advances in both theory and 
experiment for the effects of magnetism on high pressure properties. 
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Fig. 15. - Band gaps as a function of pressure for different LDA + U solutions for FeO: (a) U = 
4.6 eV, (b) U = 6.0 eV, (c) U = 8.0 eV. Solid symbols, rhombohedral solutions; open symbols, 
monoclinic solutions. Diamonds, cubic lattice; circles, rhombohedrally-strained FeO. 
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